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It has been noted a long time ago that a term of the form ^(e^/27r/i)B • E may be added to the 
standard Maxwell Lagrangian without modifying the familiar laws of electricity and magnetism. 
is known to particle physicists as the 'axion' field and whether or not it has a nonzero expectation 
value in vacuum remains a fundamental open question of the Standard Model. A key manifestation 
of the axion term is the Witten effect: a unit magnetic monopole placed inside a medium with 
^ 7^ is predicted to bind a (generally fractional) electric charge —e{0/2Ti + n) with n integer. 
Here we conduct a first test of the Witten effect based on the recently established fact that the 
axion term with — n emerges naturally in the description of the electromagnetic response of a 
new class of crystalline solids called topological insulators - materials distinguished by strong spin- 
orbit coupling and non-trivial band structures. Using a simple physical model for a topological 
insulator we demonstrate the existence of a fractional charge bound to a monopole by an explicit 
numerical calculation. We also propose a scheme for generating an 'artificial' magnetic monopole 
in a topological insulator film, that may be used to facilitate the first experimental test of Witten's 
prediction. 



I. AXIONS 



The idea of the axion was introduced in 1977 by Pec- 
cei and Quinn^ as a means to solve what is known as 
the 'strong CP problem' in the physics of strong inter- 
actions. The strong CP problem, the details of which 
are quite subtle, has to do with the vacuum structure of 
Quantum Chromodynamics. In simple physical terms it 
can be stated as a question: Why is the electric dipole 
moment of the neutron (currently unobserved) so small? 
The Standard Model predicts a value for the neutron 
dipole moment \dn\ ~ 10~^^^ e cm, with of order unity, 
that should be readily measurable. Peccei-Quinn's solu- 
tion promotes ^ to a dynamical field describing a new 
elementary particle, the axion, whose vacuum expecta- 
tion value has relaxed to a very small value, explaining 
the smallness of |dn|. The actual value of ^, and the 
validity of the Peccei-Quinn solution and its variant 
remain open questions of considerable importance to fun- 
damental physics. The axion is also believed to be a vi- 
able candidate for the elusive dark matter that comprises 
the majority of matter in our univers^ and is subject to 
active experimental searches.^^ 

In a remarkable development axion electrodynamics 
has recently emerged as a key tool in the description 
of crystalline solids called strong topological insulators 
(STIs). These three-dimensional time-reversal invariant 
(TRI) materials possess anomalous band structures char- 
acterized by a Z2- valued topological invariant.^ ^ This in- 
variant, called counts the number of topologically pro- 
tected gapless surface states (modulo 2). A non-zero in- 
variant means that the surface of such an insulator wi ll be 
metallic. This behavior has been predicted to occui^^^ 
and subsequently experimentally discovered^^^^^ in sev- 
eral 3-dimensional solids such as Bii_a^Sba; alloys and 
Bi2Se3, Bi2T e3 crystals. More recently it has been 
realized, '""^^'^^ remarkably, that the electromagnetic re- 



sponse of a STI is characterized by the axion term 

with ^ = TT, the only non-zero value permitted by the 
time-reversal symmetry. When the time-reversal symme- 
try is broken, e.g. in a crystal showing weak magnetism, 
9 can acquire an arbitrary value. Fluctuations in the 
magnetic order parameter then act as a dynamical axion 
field and can be thought of as emergent axion particles.'^ 
Thus, aside from possible practical applications, crys- 
talline solids with topologically non-trivial band struc- 
tures have the potential to provide tabletop laboratories 
for the testing and exploration of fundamental physical 
paradigms. 

A fundamental property of the axion medium is the 
Witten effect!^ in the quantum theory, a magnetic 
monopole of unit strength (i.e. projecting magnetic fiux 
<l>o = hc/e) immersed in an axion medium must carry 
electric charge —e{0/2Ti + n) with n integer. This ef- 
fect, although theoretically well established, has never 
been experimentally tested because until now both a 
suitable axion medium and the means to produce a 
magnetic monopole have been lacking. In this study 
we demonstrate how the connection between the axion 
respons^^S and strong topological insulator^^ may 
serve to overcome both obstacles. We remark that a ID 
realization of the Witten effect in antiferromagnetic spin 
chains was proposed a long time ago.^^ Here we furnish 
the first concrete physical example of the Witten effect in 
3D by modeling a STI with a magnetic monopole inserted 
in its bulk. We show that the monopole binds a fractional 
charge ±e/2 consistent with Witten's prediction. We 
then discuss possible ways to overcome the second ob- 
stacle by creating an emergent magnetic monopole in a 
topological insulator. This can be achieved by exploiting 
the degrees of freedom associated with a vortex in the 
exciton condensate that may emerge in a thin film topo- 
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logical insulator under external bias !^ We conclude that 
the prospects for experimental verification of the Wit- 
ten effect in a tabletop experiment using a STI appear 
promising. 

II. MONOPOLE AND THE WITTEN EFFECT 

We start with a brief overview of the Witten effect. 
Although the effect is quantum-mechanical in nature 
its essence can be understood by studying the classical 
Maxwell's equations modified in the presence of A>Caxion- 
The axion term revises both Gauss' law and Ampere's 
law by adding extra source terms,^^ 



V X B 



where a — /he is the fine structure constant and p, j 
are the electric charge density and current, respectively. 
We observe that for uniform, constant ^, Eqs. ([2| and 
(|3| revert to the familiar Maxwell's equations, consistent 
with the notion that A^axion can be written as a total 
derivative in this case. An important related propert}'^^ 
is the periodicity under 6 ^ 6 -\- 27rn of the axion action 
«5axion, implying that can be chosen from the interval 
[0,2^). 

Now consider a unit monopole, V • B = <l>o^(r), placed 
at the origin, in a medium initially characterized by 
^ = 0. We wish to understand what happens when we 
turn on ^ as a function of time (but keep it uniform in 
space). To this end we set = and j = (no currents 
in vacuum) and take the divergence of Eq. ^ to obtain 
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V-B = 0. 



(4) 



We see that an electric field is generated in this process. 
Integrating Eq. Q over space and time, we find that this 
field can be thought of as originating from a point electric 
charge Q located at the origin with magnitude 
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where is the net change in and we assumed 
that there was no initial electrical charge bound to the 
monopole, as should be the case for a charge-conjugation 
and parity (CP) invariant theor}W with ^ = 0. In a 
topological insulator A^ = tt, thus one expects a mag- 
netic monopole to bind fractional charge 



Q 



(6) 



The integer n accounts for the possibility of binding extra 
electrons, which can always occur - only the fractional 
part of Q is non-trivial. 



III. TOPOLOGICAL INSULATOR AS AN 
AXION MEDIUM 



We now specify our model for a topological insulator 
and show that it indeed possesses the axion term. In 
order to minimize computational difficulties we consider 
a very simple model, inspired by Ref. '22', with electrons 
hopping on the cubic lattice with two orbitals per site, 
denoted as c and d. The Hamiltonian H = Hso + Hcd, 
consists of a spin dependent part, with hopping between 
neighboring sites of the lattice. 



so 



h.c, 



(7) 



where = (cj^, Cj^, dj^^ dj^)^ ^ j labels sites of the cubic 
lattice, and are Pauli matrices in orbital and spin 
space, respectively, ji = x^y^ and the spin-independent 
terms that connect the two orbitals. 



(ij) 



h.c. 



(8) 



Although the model specified by Eqs. ^ and ([8| prob- 
ably does not describe any real solid, it is physical in 
that it is local in space and preserves time reversal and 
inversion symmetries. We show below that for a range 
of parameters (A, e, t) it represents a topological insula- 
tor, and therefore can be adiabatically deformed into any 
of the more realistic modelP'^^'^^'^ characterized by the 
same topological invariants. Any physical property that 
depends only on the topological invariants, such as the 
electrical charge bound to a monopole, can thus be cal- 
culated in the present model and the result will remain 
applicable to any topological insulator in the same topo- 
logical class. We note that a similar model was used in 
Refs. ll0ll7l25l 

Our Hamiltonian has a simple representation in mo- 
mentum space, H — ^^Hk^k with 



'Hk = -2A ^ r^cr^ sin + r^mk, 



(9) 



and mk = e — 2t cos k^. The spectrum of excitations 
has two doubly degenerate bands. 



= ±A/4A^(sin kx + sin ky + sin kz) 



(10) 



In the limit e, t ^ the bands touch at 8 non-equivalent 
Dirac points located at F^^^^^^^^^) = 7r(na^, n^, n^) with 
= 0, 1. These F^ also coincide with the 8 time-reversal 
invariant momenta^ (TRIM). When e, t are small but 
non-zero, the low-energy excitations of the system can 
be described in terms of 8 massive Dirac Hamiltonians 



(11) 



obtained by a straightforward expansion of Hi^ to lin- 
ear order in momentum in the vicinity of F^. Here 
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TABLE I: Z2 indices (z^o; 1^11^21^3) calculated according to Ref. 
[9] for our model insulator at half filling and the corresponding 
values of the axion parameter 0. is the important ^stro ng' 
invariant while z^i=i,2,3 are the so called 'weak' invariants^ 
which do not play a role in the present study but we list them 
here for completeness. It is assumed that A,t > 0. WTI 
denotes a 'weak' topological insulator. 



Parameters 


Z2 class 


Insulator type 


Axion 


|e| > 6t 


(0;000) 


trivial 





-6t < e < -2t 


(1;111) 


STI 


TT 


-2t <e < 2t 


(0;111) 


WTI 





2t <e<6t 


(1;000) 


STI 


TT 



v^^ = — 2A( — l)^'^ are the Cartesian components of the 
Dirac velocities at and = mr^ are the correspond- 
ing Dirac masses. The Hamiltonians H^. show spectra 
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The system described by Hamiltonian ^ is inversion 
symmetric and we can thus employ the method devised 
in Ref. 9 to determine the topological class of its insulat- 
ing phases when the negative-energy bands are filled with 
electrons. This straightforward method requires comput- 
ing the eigenvalues of the parity operator at the 8 TRIM 
for the occupied bands. We find four distinct phases, de- 
pending on parameters e and t, two of which are STI. 
The complete results are listed in Table I. 

According to general consider at ion4^^^ the STI 
phases should exhibit the axion term with = tt. We 
evaluate for our model using the non-abelian Berry con- 
nection A^^ = —i{a'k\di\f3'k) and the formulgP^HH 



= - f < 



(13) 



where |/3k) is an eigenstate of Hk, the trace extends over 
occupied states and di = d/dki. The integral indicated 
in Eq. ([l3| is generally difficult to evaluate and numerical 
methods must be used to obtain for an arbitrary band 
structure. For a model in which the band touching is 
described by Dirac Hamiltonians, however, a simple ana- 
lytical evaluation of Eq. (13) is possible by noticing that 
in the limit of a small Dirac mass the entire contribution 
to the integral comes from the Dirac points. We show in 
the Appendix that each Dirac point contributes 



{ill \ 
= --sgn [v^v v^rrii) 



(14) 



to the total = Y.i^ mod 27r. Although Eq. ([14| has 
been derived for the case \m£\ <C |A| we expect it to be 
more generally valid for all TRI Hamiltonians that can 
be deformed into the form of Eq. (|9| . This is because the 
value of ^ in a TRI insulator is quantized and can only 
change when a band crossing closes the gap. As long as 
band crossings occur only at points and are described 
by Dirac Hamiltonians, will be determined by Eq. (14), 
even when \m£\ is not small. 



It is easy to see that Eq. (14) gives the anticipated 
result for our model. Consider first the situation when 
e > 6t > which according to Table I is a trivial insula- 
tor. In this case all 8 masses are positive. Since four 
of the products v^VyV^ are positive and four are nega- 
tive Eq. (14) gives 6> = 0, as expected. Now suppose we 



tune e so that its value drops below 6t. This reverses the 
sign of a single Dirac mass ai i = (0,0,0). The corre- 
sponding 0^ also reverses sign and we obtain 6> = tt, as 
expected for a STI. For all other cases the values of 
are listed in Table I. These results confirm the one-to- 
one correspondence between the Z2 invariant and the 
axion parameter 6 expressed by ^ = ttuq, as expected on 
very general grounds 



IV. CHARGE BOUND TO A MONOPOLE 

We now consider a magnetic monopole in the interior 
of a STI. We model this situation by the Hamiltonian 
H defined by Eqs. ^ and ([8| with a monopole posi- 
tioned at the center of a cubic unit cell. The magnetic 
field of the monopole couples to both the electron charge 
and the electron spin through the orbital and Zeeman 
couplings, respectively. The form of the orbital coupling 
is dictated by gauge invariance and is thus universal; in 
our lattice model it is implemented by the Peierls substi- 
tution, which attaches factors e*^^^ to all hopping terms 
connecting sites i and j. Here i^ij = (27r/<l>o) J- A-d\ and 
A is the magnetic vector potential. The Zeeman coupling 
is of the form — ^/i^B ■ S/h where 113 = eh/2meC is the 
Bohr magneton and S denotes the electron spin. For free 
electrons g is close to 2 but in solids the effective g can 
be substantially larger. The Zeeman coupling thus leads 
to an additional term in the Hamiltonian, 



Hz 



(15) 



where Bj is the magnetic field at site j of the lattice. This 
term is non-universal and its importance will depend on 
the ratio of to the other relevant energy scales in 

the Hamiltonian set by A, e and t. 

We solve the Hamiltonian H = Hso + Hcd + Hz in 
a cube containing sites by exact numerical diago- 
nalization. The monopole is placed inside the central 
unit cell (at the origin), so that the magnetic field of 
the monopole is B = (<l>o/47rr^)f. We choose a gauge in 
which the system retains the four-fold rotational sym- 
metry around the z axis,^^ A = — <^>o(l + cos 0)\/(f^ with 
(6>, cp) the spherical angles. Exploiting this symmetry we 
are able to simulate system sizes up to L = 20, which 
requires diagonalizing a complex valued Hermitian ma- 
trix of size 1(4 X 20^) = 8, 000. In order to calculate the 
charge density at half filling we require knowledge of all 
the occupied eigenstates. 

We diagonalize the Hamiltonian, once with the mag- 
netic monopole and once without, obtaining charge den- 
sities pi and respectively. The monopole-induced 
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FIG. 1: (Color online) Charge density in our model TI on the cube-shaped lattice with 20"^ sites with a unit monopole at its 
center, with parameters t = A, e = 4t, leading to a bulk gap A = 4t. (a) Charge density 5p of the three closest layers below 
the monopole, for ^ = 0. (b) The excess charge 5Q{r) (in units of e) for different Zeeman coupling g. The knee feature seen at 
r = 10 corresponds to the radius at which the sphere used to calculate 5Q{r) first touches the system boundary, (c) Log-log 
plot of 5Qg — SQo showing the power-law approach ^ of the accumulated charge to its assymptotic value of 1/2. The 
least-square fit yields exponents a = 2.85, 3.04, 2.79 for = 2, 6, 10, respectively. We attribute the deviations of the numerically 
determined exponent a from the expected value of 3 to the finite size effect. 



charge density Sp = po — pi is plotted in Fig. [T^. To 
determine the total charge bound to the monopole we 
calculate the excess accumulated charge in a sphere of ra- 
dius r centered on the monopole, 5Q{r) = X]|ri|<r ^P(^i)- 
We find (Fig. [TJd) that it saturates at — e/2 to within 4 
significant digits, comparable to the accuracy of our nu- 
merics. For ^ = we find two localized zero modes, 
one at the monopole and one on the surface. Fractional 
charge bound to the monopole can be understood in this 
case by appealing to the standard arguments^^ de- 
veloped originally to describe charge fractionalization in 
polyacetylene.^^ Briefly, when a topological defect (such 
as a domain wall in polyacetylene) produces a localized 
zero mode inside the gap in a particle-hole symmetric 
system, one can show that the spectral weight of the 
state contains equal contributions from the valence and 
the conduction bands. Thus, the valence band shows a 
net deficit of half a state in the vicinity of the defect. 
This translates into the defect carrying fractional charge 
±e/2, the sign depending on whether the zero mode is 
empty or occupied. 

Like in polyacetylene we find the saturation of charge 
to be exponential ~ exp(— r/^), where ^ ~ 1/A and A 
is the bulk gap. 

When ^ > 0, the Zeeman coupling causes changes in 
the charge distribution near the monopole but the total 
accumulated charge remains quantized at —e/2. In this 
case there are no exact zero modes in the spectrum and 
SQ{r) approaches e/2 as a power law with the exponent 
close to -3 (Fig.jlj^). 

The power-law dependence can be understood as fol- 
lows. The Zeeman term acts as an additional time- 
reversal breaking field which modifies the value of ax- 
ion away from tt close to the monopole. This causes 
non- vanishing and thus, according to Eq. ([2|, addi- 
tional contribution to the effective charge density. The 
simplest assumption, SO ^ B^, gives 5p ~ • B ~ 
and SQ ~ r~^, a power law but with the exponent not 



quite in agreement with our numerical simulation. On 
further reflection one realizes that in our model 50 can- 
not be proportional to B^ but rather must be propor- 
tional to its gradients. This is because in the presence 
of a uniform Zeeman term the system remains inversion- 
symmetric. Inversion symmetry dictates quantized value 
of 6> = 0,7r even when T is explicitly broken. Thus, 
non- vanishing 60 requires spatially varying Zeeman field. 
The simplest assumption satisfying these requirements is 
\/0 ~ V^B. In the vicinity of the monopole one finds 
Sp ~ • B ~ and SQ ~ in agreement with our 
numerical results. 



We note that the above considerations are based on 
the effective axion action ([T]) and apply on lengthscales 
large compared to ^. The power law tail in the fractional 
charge distribution for ^ > appears on top of a short- 
lengthscale structure with a roughly exponential profile 
that is controlled by the properties of the microscopic 
Hamiltonian and is thus non-universal. At the interme- 
diate lengthscales the interplay of the two contributions 
can give rise to interesting structures such as the peak in 
SQ(r) at r ~ 2.5 seen in Fig. lb for g — 10. 



By the same method described above we have inves- 
tigated spin density induced by the monopole. We find 
that there is no net spin (S) attached to the monopole. 
Thus, in addition to charge fractionalization, a magnetic 
monopole inserted in a STI constitutes an example of 
spin-charge separation in three spatial dimensions. This 
is perhaps not surprising in view of the fact that spin- 
orbit coupling present in the Hamiltonian ([7|) breaks the 
SU(2) spin symmetry and, as a result, electron spin is 
not a good quantum number in the model describing our 
system. 
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V. PROPOSAL FOR EXPERIMENTAL 
REALIZATION 

Although there is no known theoretical principle 
that prohibits the existence of fundamental magnetic 
monopoles in nature}^ have been observed to date 
despite extensive searches. This null observation has 
led to a consensus that fundamental monopoles either 
do not exist for some heretofore unknown reason or they 
are very rare in our part of the universe. In either case 
the observed absence of fundamental monopoles poses a 
challenge to the idea of experimental verification of the 
Witten effect using a STL At best, one could conceive of a 
new scheme for possible detection of magnetic monopoles 
using a STI in the role of a sensor if a convenient way to 
detect the fractional charge could be found. 

A much more promising avenue for the verification of 
the Witten effect is suggested by exploiting emergent 
instead of fundamental monopoles. A classic example 
of such an emergent behavior in a crystalline solid is 
the 2007 theoretical predictioii^ and the subsequent ex- 
perimental observation^^ of monopoles in frustrated 
magnetic systems called 'spin ice', realized in certain 
magnetic pyrochlore compounds such as Dy2Ti207 or 
Ho2Ti207. Magnetic monopoles in these systems arise as 
elementary excitations above the collective ground state 
of spins and the monopole-like magnetic field configura- 
tion originates from the magnetic moments of the con- 
stituent spins. In principle, the emergent monopoles in 
the spin ice could be used to test the Witten effect if a 
compound that is simultaneously a STI and a spin ice 
could be identified. Unfortunately no such material is 
known at present although we note that STI behavior 
has been theoretically predicted to occur in crystals with 
the same pyrochlore structure^^ that underlies the spin 
ice behavior. It is thus possible that a suitable material 
will be discovered in the future. 

Here we focus on a different type of emergent magnetic 
monopole that can arise in a thin film STI placed in a 
uniform external electric field. The basic idea and the 
feasibility of its experimental realization have been dis- 
cussed in Ref. 20, Following that work we envision the 
simplest STI with just one gapless Dirac state per surface 
and the chemical potential initially tuned to the neutral 
point. When a strong enough electric field is applied per- 
pendicular to the plane of the film the chemical poten- 
tial undergoes a shift that is opposite in the two surfaces. 
This creates a small electron Fermi surface in one surface 
and a small hole Fermi surface in the other. The essence 
of the proposal^^ lies in the observation that an arbitrar- 
ily weak Coulomb interaction between the surface states 
produces an exciton condensate, which may be viewed 
as a coherent fiuid of electron-hole pairs drawn from the 
opposite surfaces. Such an exciton condensate is char- 
acterized by a complex scalar order parameter which 
can fiuctuate in space and time. In particular ^ = ^oe*^ 
can contain vortices - point-like topological defects with 
the phase x winding by ±27r around a vortex. It has been 



pointed out in Ref. |20| that to electrons in a STI such a 
vortex is indistinguishable from a 'planar monopole', i.e. 
a monopole with magnetic field radiating in the plane of 
the surface. 

A planar monopole can be viewed as an adiabatic de- 
formation of an ordinary monopole achieved by fiattening 
the field lines in a cylindrically symmetric fashion. One 
expects that the total charge bound to the monopole via 
the Witten effect should be insensitive to such an adia- 
batic deformation and therefore a vortex in the exciton 
condensate should bind fractional charge ±e/2. This in- 
deed has been argued to happen in Ref. 20 based on the 
Dirac equation describing the low-energy physics of the 
surface states in the presence of the exciton condensate. 
Here, taking a more general point of view, we establish 
the existence of the fractional charge in such a condensate 
by studying a planar monopole embedded inside a STI. 
Our calculation below does not rely on the low-energy 
approximation for the surface states and is insensitive 
to the detailed microscopic structure of the condensate. 
Rather, it exploits only the most fundamental property 
of the STI given by its nontrivial axion response. 

In general, the fractional charge is expected to be ro- 
bust against weak disorder that does not break TRI. Such 
a disorder will be present in a real sample and we model 
it here by adding a term 

hd = J2 + E -"J^i^-^i' (16) 

j 3 

to the Hamiltonian. The first term represents a parity- 
preserving on-site disorder (independent of the orbital), 
while the second term is a parity breaking disorder. 

As before we solve the Hamiltonian 
H = Hso + Hcd + H]j^ in a cube containing sites, by 
exact numerical diagonalization. The planar monopole 
projects an effective magnetic field Beff = {^o/27rr)d{z)r 
(in cylindrical coordinates) and the vector potential 
can be chosen as A = {^o/27T)(p5{z)z. The effective 
field does not couple to electron spin,!^ so there is no 
Zeeman term in this case. The disorder coefficients /i^ 
are chosen from a Gaussian distribution with standard 
deviation /i. Note that the disorder breaks the four-fold 
rotational symmetry of the system, so we cannot exploit 
this symmetry in this case to efficiently diagonalize the 
Hamiltonian. Consequently we are limited to system 
sizes up to L = 14. For weak disorder /i <C A the charge 
bound to the planar monopole remains — e/2 (see Fig.[|), 
and for strong disorder /j. A the charge bound is zero. 
Remarkably, even for fairly significant disorder (such 
that it generates charge density fiuctuations comparable 
to the charge density induced by the monopole) the 
difference in charge density Sp shown in Fig. [2]P is only 
weakly affected. 

In the framework of the current proposal the key ingre- 
dient required to produce a monopole-like configuration 
is the exciton condensate. As explained in Ref. 20 it 
is difficult to reliably estimate the critical temperature 
Tec for the formation of the exciton condensate, but un- 
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FIG. 2: (Color online) A cubic sample of a TI including disorder with a planar unit monopole at its center, size L = 14 and 
parameters as in Fig. [1] In all cases shown we use the same disorder realization but vary its overall strength parametrized by 
/i. Panels (a) and (b) show charge density 2 — pi for the layer just below the planar monopole for weak disorder /i = 0.05A, 
and (c) larger disorder /x = 0.20A. Panels (e) and (f) show the difference in charge density 5p — pQ— pi for /x/A = 0.05,0.20, 
respectively, for the same layer, (d) The excess charge 5Q(r) (in units of e) for different values of disorder strength p. The 
inset shows a close up of the saturation. At this scale a small deviation from the expected asymptotic value 1/2 that increases 
with the disorder strength becomes visible. We attribute this deviation to the finite-size effect in our numerical calculation. 
This identification is supported by the fact that the deviations grow more pronounced for smaller system sizes and close to the 
surface. Also, it is consistent with the notion that the bound charge is localized on the length scale ^ ^ 1/A which increases 
as the disorder reduces the spectral gap. 



der optimal conditions it should be higher than it is in 
bilayer graphene, where the occurrence of this effect is 
hotly debated. Once the exciton condensate is formed, 
vortices can be nucleated by applying an in-plane mag- 
netic field. Since the exciton condensate is itself insulat- 
ing, the main conduction channel in this situation will 
be through vortices, each carrying — e/2 charge. Frac- 
tional charge of the carriers then can be detected using 
established techniques. -^^-^ 



VI. OUTLOOK AND OPEN QUESTIONS 

Predicted more than 30 years ago in the context of 
high-energy physics, but never before observed in a real 
or numerical experiment, the Witten effect is realized in a 
strong topological insulator. A unit magnetic monopole 
inserted in a model STI binds electric charge —e/2 in 
accordance with the prediction^^ and furnishes a rare ex- 
ample of charge fractionalization and spin-charge sepa- 
ration in 3 spatial dimensions. In the special case when 
the underlying system possesses particle-hole symmetry 
and when the Zeeman term ( 15 ) can be neglected, the ap- 



pearance of fractional charge follows from the same 'zero- 
mode' arguments that underlie charge fractionalization in 
a one-dimensional system of fermions coupled to a scalar 



field with a soliton profil#^^ as realized in dimerized 
poly acetylene.^ In the more general case when the Zee- 
man term or weak disorder are present, there exist no ex- 
act zero modes in the spectrum of electrons yet the frac- 
tional charge remains precisely quantized. This reflects 
the more subtle topological order that underlies the axion 
response of a STI, which is robust against any weak per - 
turbation that respects time reversal symmetry.l^EE^EIl 

An interesting open question is the fate of the Witten 
effect in the presence of magnetic disorder. Experimen- 
tally, this can be implemented by adding a small con- 
centration of magnetic ions (such as Fe or Mn) into a 
topological insulator. Although at the microscopic level 
9 is no longer quantized in the presence of T-breaking 
perturbations, we expect its effective value, relevant to 
the physics at long lengthscales, to remain pinned at tt 
as long as the magnetic moments stay disordered. This 
is because the net magnetic moment in a macroscopic 
region containing many impurities will effectively vanish. 
Qualitatively, this suggests that the Witten effect may 
survive inclusion of a moderate concentration of mag- 
netic dopants with randomly oriented moments. At low 
temperatures moments may order ferromagnetically.^In 
this case both T and the inversion symmetry are bro- 
ken (the latter due to the random position of magnetic 
dopants) and the effective can acquire an arbitrary 
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value. In this situation we expect a monopole to still 
bind fractional charge according to Eq. ([6| but we leave 
a detailed study of this case to future investigation. 

Can the Witten effect be observed experimentally in 
the near future? We believe that the answer is affir- 
mative. One essential ingredient, the axion medium, is 
now widely available in any of the recently discovered 
STIs.'^^K^ If an emergent monopole can be realized, ex- 
ploiting the propo sed exciton condensate, the spin ice- 
type physics,^^^^^ or by some other means, then the ex- 
perimental challenge is reduced to designing a suitable 
method for the detection of fractional charge bound to 
the monopole. The fractional charge of elementary exci- 
tations in fractional quantum Hall fluids has been previ- 
ously detectecP^l^ and it should be possible to adapt 
these methods to topological insulators. In this way, 
studies of crystalline quantum matter with non-trivial 
topological properties could help settle one of the endur- 
ing challenges of fundamental physics and provide new 
insights into the behavior of electrons placed in unusual 
situations. 
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Appendix A: Evaluation of d 

The feature that makes our model Hamiltonian ([9| 
easy to analyze, is its matrix structure in the combined 
orbital/spin space, which consists of four anticommut- 
ing 4x4 matrices. Such matrices are by convention de- 
noted 7^, \i = 0, 1, 2, 3 and form a representation of the 
Clifford algebra defined by the anticommutation relation 
{7/i,7i/} = 2(5^i,. In terms of these matrices Eq. ([9| can 
be written as 



'^k = ^7/i^/i(k) 



(Al) 



with (i^(k) = {rn\^^ —2X sin kx, —2X sin ky, —2X sin kz). 
The actual form of our 7 matrices is apparent by com- 
paring Eq.(Al) to (|9| but our result for is independent 
of any particular representation, as long as the matrices 
satisfy the requisite commutation relation. In fact in the 
subsequent calculations it will be advantageous to use a 



nii = DodiDi 
n2i = DsdiDi 
nsi = D2diDi 



DidiDo ^ DsdiD2 - D2diDs, 
DidiDs + D2^^Do - DodiD2, 
DidiD2 + Do^^Ds - DsdiDo, 



with Djj^ = d^/y/2E. After substituting these into Eq. 
(13) a tedious but ultimately straightforward calculation 



leads to the expression 

27r Jbz E 

Using the values of d^ for our model given below Eq. 
(Al), one obtains a complicated integrand in terms of 
trigonometric functions. Although the value of the inte- 
gral is guaranteed to be either or tt in practice it is not 
obvious how to perform the required three-dimensional 
integration. However, it is clear from the structure of 



the integrand in (A3) that in the limit |mk| <C A contri- 



butions to the integral come only from the vicinity of the 
8 Dirac points. We evaluate these 8 contributions sepa- 
rately by linearizing d's as <i^(k) = (m^, v^kx^Vyky, ^i^^z) 
and obtain 

27r J {AX^k^ + m|)2 

An elementary evaluation then yields the result quoted 
in Eq. (pi). 
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